Theory of ground state factorization in quantum cooperative systems 
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We introduce a general analytic approach to the study of factorization points and factorized ground states 
in quantum cooperative systems. The method allows to determine rigorously existence, location, and exact 
form of separable ground states in a large variety of, generally non-exactly solvable, spin models belonging 
to different universality classes. The theory applies to translationally invariant systems, irrespective of spatial 
dimensionality, and for spin-spin interactions of arbitrary range. 
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Introduction- Quantum engineering and quantum techno- 
logy have been developing at a fast pace in recent years. 
Quantum devices are being vigorously pursued for applica- 
tions ranging from nano-sciences to^juantum information and 
entanglement-enhanced metrology [1]. Despite a large variety 
of possible implementations involving different physical sys- 
tems, many relevant properties of such devices can be investi- 
gated in a unified setting by appropriate mappings to quantum 
spin models 12j,|3[]. Thus, control of ground state entanglement 
in quantum spin systemsplays an important role in quantum 
technology applications 14J]. On the other hand, knowledge of 
exact solutions endowed with precisely determined properties 
of separability or entanglement, can be of great relevance in 
the study of advanced models of condensed matter and coop- 
erative systems that are in general not exactly solvable. 

The occurrence of totally factorized (unentangled) ground 
states of quantum many-body systems was first discovered 
in the one-dimensional anisotropic Heisenberg model with 
nearest-neighbor interactions |5[]. This result was later re- 
derived and extended to two dimensions using quantum 
Monte Carlo numerical methods jg]. Complex quantum sys- 
tems exhibiting cooperative behaviors, whose ground states 
are typically entangled fl7[], may thus admit, for some non triv- 
ial values of the Hamiltonian parameters, a ground state which 
is completely separable. The phenomenon of ground state 
factorization appears to be associated with the presence of 
an "entanglement phase transition" with no classical counter- 
part |8|] ; furthermore, for the purposes of quantum engineering 
applications that employ distributed entanglement in order to 
manipulate and transfer information ||9[], factorization points 
need to be exactly identified and avoided to guarantee the reli- 
able implementation of quantum devices. Finally, for models 
not admitting exact general solutions, achieving knowledge of 
the exact ground state, even if only for the restricted nontrivial 
set of parameters associated to factorization, would allow ( i) 
to prove the existence of an ordered phase and characterize it; 
(ii) to build variational or perturbative approximations around 
the exact factorized solution, that may then be used as test 
benchmarks for the validity and the precision of numerical al- 
gorithms and simulations. Unfortunately, to date, it has been 



extremely hard to go beyond the pioneering result of Kurmann 
et al. ]5[] despite the fact that, to prove total factorization, it 
would suffice to show the vanishing of the von Neumann en- 
tropy of entanglement, or of the linear entropy (tangle) II 1 Oil - 
The difficulty resides in the fact that, with few special excep- 
tions, explicit analytic expressions for these measures of en- 
tanglement cannot be obtained. Hence, the only possibility 
to gain insight on factorization in systems of increasing com- 
plexity has relied so far on heuristic or numerical approaches. 

In the present work we introduce a general analytic method 
that allows to determine exactly the existence of factorized 
ground states and to characterize their properties in quantum 
spin models defined on regular lattices, in any spatial dimen- 
sion, and with spin-spin interactions of arbitrary range. In 
correspondence to rigorously established ground state factor- 
izability, the method also allows to determine novel sets of 
exact solutions in generally non exactly solvable models. Pre- 
vious particular findings for models with short range interac- 
tions are rigorously re-derived and extended within a unified 
framework inspired by concepts of quantum information sci- 
ence. The method is built on a formalism of single-spin, or 
single-qubit, unitary operations (SQUOs) and associated en- 
tanglement excitation energies (EXEs), previously introduced 
for the characterization and quantification of entanglement in 
systems of quantum information llll Il2ll . The novel tech- 
niques exploit the fundamental property enjoyed by the EXEs, 
of vanishing if and only if a pure state is fully factorized 112 ] . 
For any given Hamiltonian, the strategy to the understanding 
of factorization is first to assume as working point a phase 
endowed with some kind of magnetic order. Next, by impos- 
ing the vanishing of the EXE and of the linear entropy, one 
derives a closed set of conditions whose solutions determine 
uniquely the occurrence (or the non occurrence) of factoriza- 
tion points at which a quantum ground state is completely dis- 
entangled. Besides the rigorous determination of novel fac- 
torization points and exact solutions of generic quantum spin 
models, the method allows as well to re-derive analytically 
the fewpreviously known results on ground state factoriza- 
tion JUSH. 



The method- To fix ideas and notations, let us consider 
general, translationally invariant, exchange Hamiltonians for 
spin- 1/2 systems on d-dimensional regular lattices, with spin- 
spin interactions of arbitrary range and arbitrary anisotropic 
couplings. This class of Hamiltonians encompasses a very 
large set of models describing different spin systems and span- 
ning several universality classes like, among others, the Ising, 
XY, Heisenberg, and XYZ symmetries. The general Hamilto- 
nian can be written in the form 
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Here i (and similarly I) is a d-dimensional index vector iden- 
tifying a site in the lattice, Sf (a = x,y, z) stands for the 
spin- 1/2 operator on site i, h is external field directed along 
the z direction, r = \i — l\ is the distance between two lat- 
tice sites, and J r a is the spin-spin coupling along the a direc- 
tion; translational invariance implies that it depends only on 
the distance r between the spins. Without loss of generality, 
one can impose \J X \ > | JI|, | J z \ W. This condition guaran- 
tees that at a particular value of the external field h = h c the 
system undergoes a quantum phase transition at zero temper- 
ature: For h < h c the system is in an ordered phase which 
may, or may not, be associated to a non-vanishing order pa- 
rameter corresponding to the ground-state expectation of S%. 
The order parameter M x = (S%) in the case of ferromagnetic 
order, and M x = (— 1)'(5^) in the anti-ferromagnetic case. 
In the following, we specialize to the anti-ferromagnetic case; 
trivial modifications are needed in the ferromagnetic case. 

Single-Qubit Unitary Operations (SQUOs) are unitary 
transformations Uk that leave all spins unaffected but for an 
arbitrarily chosen one, say, at site k, on which the SQUOs 
act as unitary, Hermitian, and traceless operators ll III . One 
can prove that there exists an element of this class, the Ex- 
tremal SQUO (E-SQUO) Uk, such that the squared Euclidean 
distance between a state \^>) and its image U$&) under the 
action of the E-SQUO coincides with the linear entropy (tan- 
gle) r = 2(1 — Tr[p1)), where pk denotes the reduced state 
of spin k lid lllll . This entanglement monotone quantifies 
the entanglement existing in state \ty) between the single-spin 
block k and the remainder of the system. If | \P) is the ground 
state of a quantum mechanical Hamiltonian H, the E-SQUO 
is uniquely associated to the aforementioned EXE, defined as 
AE = (UkHUk) -(H). A crucial property of the EXE is 
that if H is translationally invariant and [H, Uk} ^ for all 
SQUOs Uk, then the ground state is completely factorized 
if and only if AE = II 1211 . The generic spin- 1/2 models 
Eq. ([T), as well as many others, satisfy this condition. The 
above theorem can then be applied to identify the occurrence 
of factorization points. By definition, the E-SQUO can be 



written as Uk = &)i+k -*-! ® ®k. wnere ' following Ref. lllll . 
Ok can be written as 



Of =S%cose±S%sint 



(2) 



where the ± sign discriminates the two sublattices, corre- 
sponding to the sign of the staggered magnetization (SJ: ). 

Let us assume that H, Eq. (HJ, admits a factorized ground 
state. Applying the E-SQUO, Eq. (01, to the ground state, ex- 
ploiting the fact that for fully factorized states all correlations 
separate in products of single-site expectations, and imposing 
the condition AE = 0, one has that factorization requires the 
simultaneous occurrence of 



tan0 = 



Mx 
M z 



tan 6 



J X M X 
J Z M Z - h f 



(3) 



where hf is the factorizing field, i.e. the value of the external 
field for which ground state factorization occurs. The quanti- 
ties J a are the net interactions reflecting the type of magnetic 
order that exists along different axes. In the antiferromagnetic 
case, J x = J2T=i(~^y z rJ r x and J z = Y^°=i z rJ r z , where 
Z r denotes the number of sites placed at distance r from a 
given spin. Conditions (fJJ and the vanishing of the tangle 
t = 1 — 4(Mj + M z ) yield a closed expression for the phase 
8 as a function of the Hamiltonian parameters and of hf. 
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Eq. (HJi determines, independently of the magnetizations, the 
form of the candidate factorized ground state \^f) Il4ll : 



I*/ 
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(5) 



where \ipf ) are the eigenvectors of O^ with eigenvalue 1/2. 
This far, we have determined the general expression that 
a factorized ground state must assume. We are left to estab- 
lish the conditions for its existence, i.e. the conditions under 
which a state of the form Eq. (0 is indeed the eigenstate of H, 
Eq. (HJ, with the lowest energy. For each pair of spins i and j 
we introduce the pair Hamiltonian 



H. 



LI 
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where h r f is defined by the relation 1h r r = cos#(JJ — 
(—1 )'"</£). It is immediate to verify that by re-summing the 
operators Hjj over all spin pairs, one reobtains Eq. ([T), with 
h = hf. Hence, proving that \^>f) is a simultaneous eigen- 
state of all pair Hamiltonians H^, implies that it is an eigen- 
state of the total Hamiltonian H as well. Inserting the expres- 
sion of \^ff) in Eq. © yields a set of conditions that must 
be satisfied to ensure that \^>f) is an eigenstate of every pair 



Hamiltonian H,, 



LI- 



-J , y +cos 2 9J , x + (-l) r sm 2 0J , z ^O Vr, (7) 

where 9 is given by Eq. ©. By summing over r, term by term, 
all the relations in Eq. (0, and solving for hf, we eventually 
obtain the exact, general expression of the factorizing field as 
a function of the net interactions along the different axes: 



2\ {Jx 



Uz) \*Jy <J z 



(8) 



In Eq. dH) the net interaction along the y axis depends on the 
anti-ferromagnetic order on the x axis and, hence, it is given 
by J y = ^2{—l) r Z r Jy. We remark that Eq. ([8} is completely 
general and holds for lattices of arbitrary spatial dimension 
and for interactions of arbitrary range. 

We are left to determine the conditions under which the 
factorized eigenstate \*f?f) is associated to the lowest energy 
eigenvalue. A general sufficient condition is that every two- 
site reduced state, obtained from \*Bf) by a partial trace over 
all sites except the pair {i, j} (which is still a pure state since 
\^f) is factorized), is the ground state of Hjj, for every pair 
{i, j}. Given the relation existing between all the pair Hamil- 
tonians Hjj and the total Hamiltonian H, it follows that if a 
state is associated to the lowest eigenvalue of every Hjj, then 
it is associated to the lowest eigenvalue of H. To proceed, 
we need to distinguish between the various possible cases, de- 
pending on the structure of the couplings J T a . 

Models with short range interactions- By short range, or 
nearest-neighbor, we mean J r a = for all a and r > 2. 
Anti-ferromagnetic order is ensured by having J* = 1 and 
\Jy z \ < 1. For one-dimensional models, we have Z\ = 2, 
J x = —2, J y = —2Jy, and J z = 2J\. Inserting these quan- 
tities in Eqs. ([8]|U[5]l we obtain the expressions for and for 
the factorized ground state |\&/). The explicit expression of 

the factorizing field is h { f =1) = J(1 + J})(J^ ~+j[) and 

the energy per site reads g( d=1 ) = (1/8)(J Z - J y + 2) = 
(1/4) (1 + J\ + Jy). Besides reproducing the original results 
of Kurmann et al. 05D, the analytic method allows to establish 
that ground state factorization occurs for a much larger range 
of values of the couplings. Our general framework singles out 
novel instances of classical-like ground states already in this 
simple model, as pictorially sketched in Fig.[T] 




/, 



FIG. 1: (Color online) Spin-1/2 models with short range interac- 
tions: Domains in the space of couplings. Region I (vertical lines): 
Domain of parameters for which ground state factorization was origi- 
nally identified |5fl, and is established rigorously in the present paper. 
Region II (horizontal lines): Domain in which factorization is estab- 
lished for the first time, and rigorously, in the present paper. Region 
III (no lines): Domain in which factorization is rigorously excluded. 
All plotted quantities are dimensionless. 

The analytic approach can be extended immediately to 
higher-dimensional lattices. For a square lattice, the number 
of nearest neighbors is Z\ = 4, and hence the net interactions 



read J x = -4, J y = -U\, and J z = U\. Therefore, 
moving from one to two dimensions, we find that factoriza- 
tion occurs in the same domains of couplings and for the same 
value of 6, but at a factorizing field and with an energy per site 
that are twice the corresponding quantities in one dimension: 



2/i ( , d_1) ande (d=2) = 2e^ d=1) . These exact results 



■ (d=2) 

,i f - m> s 

recover, confirm, and extend recent numerical findings 16fl. 

As the resources for numerical simulations scale with the 
lattice dimension, and in the absence of analytic approaches, 
it is not surprising that no study of factorization in three- 
dimensional systems was attempted so far. However, exploit- 
ing our novel analytic method, such a study can be carried out 
exactly and straightforwardly. Namely, moving from square 
to cubic lattices we only need to insert the correct value of 
the coordination number Z\ = 6. It is then straightforward 
to prove that for three-dimensional models ground state fac- 
torization occurs at the same values of the couplings as in one 
dimension, but at a factorization field and with an energy per 
site that are three times the corresponding quantities in one 



dimension: h 



(d=3) 

/ 



34 d=1) and £ ( d = 3 )=3 £ ( d = 1 ). 



The ferromagnetic counterparts of the anti-ferromagnetic 
models can be immediately recovered performing ir/2- 
rotation around the z axis at each lattice site in either one 
of the two sublattices. Hence, given a set of couplings 
(1, J}. , J\ ) for which an anti-ferromagnetic factorized ground 
state occurs, there exists a corresponding set of couplings 
(— 1, — Jy, J*) for which a ferromagnetic factorized ground 
state occurs at the same value of the factorizing field hf, the 
same energy per site, and the same value of 9. 

Models with finite range interactions- In the case of spin 
systems with short range interactions, the set of equations 
Eq. (0 is automatically verified by any model that admits a 
real solution for Eq. (01 and Eq. ^ IU5I1 . This redundancy is 
removed as soon as one considers interactions of longer spa- 
tial range, because in these cases the net interactions do not 
depend on the nearest-neighbor couplings alone. To illustrate 
this important point, let us consider ferromagnetic models in 
arbitrary spatial dimensions and with non-vanishing interac- 
tions up to a certain finite distance s: J£ < V?' < s, 
Jx y = ^ r — s ' an d Jl = ^ r - ft i s straightforward 
to verify that a sufficient condition for ground state factor- 
ization is that the ratio of the non vanishing couplings must 
satisfy the relation J£/ ' J\ = Jy/Jy = 7r Vr < s. Otherwise, 
Eqs. (|7]i do not admit solutions. If this condition is satisfied, 
we have that the net interactions (that in the ferromagnetic in- 
stance read J x . y = J2 r Z ^l, y ) aie Jx, v = ^ J x,y where 
r = Y^,r=i Zrlr- Ground state factorization occurs at hf = 
T{JlJlf/ 2 , with an energy per site e = [T{J^. + J,J)]/8. 
Therefore, at exactly defined ratios of the couplings, systems 
with finite range interactions admit fully separable ground 
states analogous to the ones arising in the case of systems with 
only nearest-neighbor couplings. The factorizing field and the 
energy per site are increased exactly by a factor T/2 with re- 
spect to the case of models with short range interactions. 



Models with infinite range interactions- A very interesting 
limiting case is given by models with infinite range interac- 
tions, such as the fully connected or Lipkin-Meshkov-Glick 
(LMG) model 11611 . which is obtained in the limit of diverg- 
ing s and couplings of the form J r x = J xy = 2A x/y /N Vr. 
The scaling with the number of lattice sites N ensures that 
the net interactions converge to a finite value in the thermody- 
namic limit: J XtV — > 2A XtV . It is then rather straightforward 
to solve Eqs. © exactly and prove rigorously that the ground 
state of the LMG model is a fully factorized ferromagnetic 
state for A x = 1; < A y < 1; hf = \fA v at 9 = arccos hf\ 
and an energy per site e = (1 + A,,") I A. These results justify 
rigorously recent numerical findings 11311 . 

Comments.- We have introduced a simple and general ana- 
lytic approach to the exact determination of factorized ground 
states in quantum spin systems. We have applied the scheme 
to spin-1/2 models with general anisotropic Heisenberg-like 
interactions of arbitrary range and for lattices of arbitrary di- 
mensions. Besides the rigorous derivation of the few known, 
mainly numerical, results, we have showed that our method al- 
lows to determine exactly novel classes of factorization points 
in various models, generally non-exactly solvable, for dif- 
ferent lattice dimensions and for different interaction ranges. 
These novel exact solutions of non exactly solvable models 
are obtained for nontrivial sets of values of the Hamiltonian 
parameters. Furthermore, according to the general theorem by 
Kurmann et al. on factorization |5[], given any Hamiltonian 
of the form Eq. (Q]) with generic spin S > 1/2, the ground 
state of the system is factorized at the same value of the ex- 
ternal field h — hf [Eq. dH)], at which factorization occurs 
in the corresponding spin-1/2 model. Therefore, the method 
and the results derived in the present paper are straightfor- 
wardly generalized to interacting systems with arbitrary value 
of the spin which are endowed with the same Hamiltonian 
structures as in the spin-1/2 case. Further applications to 
other systems, defined, e.g., on ladders and coupled planes, 
or to models with frustration and in complex geometries, can 
be in principle carried out by suitably adapting and special- 
izing the general framework introduced in the present work. 
From a conceptual standpoint, the method realizes a rigorous 
and analytic implementation of concepts motivated by quan- 
tum information theory to obtain genuinely new insights on 
founding open questions of condensed matter physics. 
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